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We use perturbative Symanzik improvement to create a new staggered-quark action (HISQ) that 
has greatly reduced one-loop taste-exchange errors, no tree-level order a 2 errors, and no tree-level 
order (am) 4 errors to leading order in the quark's velocity v/c. We demonstrate with simulations 
that the resulting action has taste-exchange interactions that are at least 3-4 times smaller than the 
widely used ASQTAD action. We show how to estimate errors due to taste exchange by comparing 
ASQTAD and HISQ simulations, and demonstrate with simulations that such errors are no more 
than 1% when HISQ is used for light quarks at lattice spacings of 1/10 fm or less. The suppression 
of (am) 4 errors also makes HISQ the most accurate discretization currently available for simulating 
c quarks. We demonstrate this in a new analysis of the tp — n c mass splitting using the HISQ 
action on lattices where am c = 0.43 and 0.66, with full-QCD gluon configurations (from MILC). 
We obtain a result of 111(5) MeV which compares well with experiment. We discuss applications of 
this formalism to D physics and present our first high- precision results for D s mesons. 

PACS numbers: 11.15.Ha,12.38.Aw,12.38.Gc 



I. INTRODUCTION 

The reintroduction of the staggered-quark discretiza- 
tion in recent years has transformed lattice quantum 
chromodynamics, making accurate calculations of a wide 
variety of important nonperturbative quantities possible 
for the first time in the history of the strong interac- 
tion P, 0, OJ 0, HJ. Staggered quarks were introduced 
thirty years ago, but unusually large discretization errors, 
proportional to a 2 where a is the lattice spacing, made 
them useless for accurate simulations. It was only in the 
late 1990s that we discovered how to remove the lead- 
ing errors, and the result was one of the most accurate 
discretizations in use today. Staggered quarks are faster 
to simulate than other discretizations, and as a result 
have allowed us for the first time to incorporate (nearly) 
realistic light-quark vacuum polarization into our simula- 
tions. This makes high-precision simulations, with errors 
of order a few percent, possible for the first time (see jlj 
for a more detailed discussion). In this paper we present 
a new discretization that is substantially more accurate 
than the improved discretization currently in use. With 
this new formalism accurate simulations will be possible 
at even larger lattice spacings, further reducing simula- 
tion costs. 

The 0(a 2 ) discretization errors in staggered quarks 
have two sources. One is the usual error associated with 
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discretizing the derivatives in the quark action. The 
correction for this error is standard and was known in 
the 1980s |fj. The second source was missed for al- 
most a decade. It results from an unusual property 
of the staggered-quark discretization: the lattice quark 
field creates four identical flavors or tastes of quark 
rather than one. (We refer to these unphysical flavors 
as "tastes" to avoid confusion with the usual quark fla- 
vors, which are not identical since quarks of different fla- 
vor have different masses.) The missing a 2 error was as- 
sociated with taste-exchange interactions, where taste is 
transferred from one quark line to another in quark-quark 
scattering. The generic correction for this type of error 
involves adding four-quark operators to the discretized 
action. 

Taste is unphysical. Instead of one 7r + , for exam- 
ple, one has sixteen. Taste is easily removed in simu- 
lations provided that different tastes are exactly equiva- 
lent, which they are, and provided there are no interac- 
tions that mix hadrons of different taste. Taste-exchange 
interactions, however, cause such mixing, and therefore it 
is crucial that we understand such interactions and sup- 
press them. Furthermore, past experience suggests that 
the errors due to residual taste exchange are the largest 
remaining a 2 errors in current simulations. 

The corrections to the lattice action that suppress taste 
exchange were missed initially because four-quark op- 
erators are not usually needed to correct discretization 
errors in lowest order perturbation theory (that is, at 
tree level). They were first discovered empirically from 
simulations in which the gluon fields in the quark ac- 
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tion were replaced by smeared fields, which happen to 
suppress taste-changing quark-quark scattering ampli- 
tudes • This discovery led to a proper understanding 
of the taste- exchange mechanism 0, [tj , and to the first 
correct analysis of the a 2 errors in staggered-quark for- 
malisms H . The resulting "ASQTAD" quark action 
has provided the basis for almost all simulations to date 
that include (nearly) realistic light-quark vacuum polar- 
ization. 

While it is possible to remove all tree-level taste 
exchange by smearing the gluon fields appropriately, 
higher-order corrections can only be removed by adding 
four-quark operators. Recent numerical experiments sug- 
gest, however, that even these corrections can be sig- 
nificantly suppressed through additional smearing . 
These studies, while very important, were limited in two 
ways. First, they used the mass splittings between the 
sixteen tastes of pion as a probe for taste exchange; there 
are many taste-exchange interactions that do not affect 
this spectrum. Second, smearing the gluon fields intro- 
duces other types of 0(a 2 ) error which undo the advan- 
tage of an improved discretization (that is, accuracy at 
large values of a). 

In this paper we present the first rigorous analysis of 
one-loop taste exchange. Based upon this analysis we 
develop a simple smearing scheme that suppresses one- 
loop corrections by an order of magnitude on average. 
Our smearing scheme is simpler than that in ^lj and, 
unlike that scheme, it does not introduce new a 2 errors. 
The result is a new discretization for highly-improved 
staggered quarks which we refer to by the acronym HISQ. 

An added advantage of a highly improved action is that 
it can be used to simulate c quarks. Heavy quarks are dif- 
ficult to simulate using standard discretizations because 
discretization errors are large unless am <C 1, where a is 
the lattice spacing and m the quark mass. To simulate 
b quarks, for example, one would require lattice spacings 
substantially smaller than 1/20 fm — much too small to 
be practical today. Consequently high-quality simula- 
tions of b quarks rely upon rigorously defined effective 
field theories, like nonrelativistic QCD > that remove 
the rest mass from the quark's energy. While this ap- 
proach works very well for b quarks, it is less successful 
for c quarks because the quark's mass is much smaller, 
and consequently c quarks are much less nonrelativistic. 
The smaller quark mass, however, means that am c w 1/2 
for lattice spacings of order 1/10 fm, which are typical of 
simulations today. As we demonstrate in this paper, we 
can obtain few-percent accuracy for charm quarks when 
am c sa 1/2 provided we use a highly improved relativistic 
discretization like HISQ. 

In Section II we review the (perturbative) origins 
of taste- exchange interactions in staggered quarks, and 
their removal at lowest order (tree level) in perturbation 
theory, using the Symanzik improvement procedure. We 
also extend the traditional Symanzik analysis to cancel 
all (am) 4 errors to leading order in the quark's velocity, 
v/c (in units of the speed of light c). The additional 



correction needed is negligible for light quarks, but sig- 
nificantly enhances the precision of c-quark simulations. 

In Section III, we extend our analysis of taste exchange 
to one-loop order, and show how to suppress such contri- 
butions by approximately an order of magnitude, thereby 
effectively removing one-loop taste exchange from the 
theory. The resulting quark action (HISQ) is the first 
staggered-quark discretization that is both free of lattice 
artifacts through 0(a 2 ) and effectively free of all taste 
exchange through one-loop order. We also identify the 
order a s (am) 2 corrections to the action that are impor- 
tant for high-precision c physics. 

We illustrate the utility of our new formalism, in Sec- 
tion IV, first by examining its effect on taste splittings in 
the pion spectrum. We also show how to directly measure 
the size of taste-exchange interactions in simulations, by 
comparing ASQTAD and HISQ simulations. There has 
been much discussion recently about the formal signifi- 
cance of taste-exchange interactions With our pro- 
cedure it is possible to directly measure the effect of these 
interactions on physical quantities. We demonstrate the 
procedure by showing that residual taste-exchange inter- 
actions in HISQ contribute less than 1% at current lattice 
spacings to light-quark quantities such as meson masses 
and meson decay constants. 

We report on a new simulation of c quarks and charmo- 
nium using the HISQ formalism in Section V. This is the 
most severe test of the HISQ formalism since am c w 1/2 
for the lattice spacings we use, but we show that the for- 
malism delivers results that are accurate to a couple of 
percent, making it the most accurate formalism there is 
for simulating c quarks. We demonstrate this in a new 
high-precision determination of the tp—rjc splitting. This 
analysis confirms our expectation that the largest a 2 er- 
rors in the ASQTAD action are associated with taste 
exchange, and we show that these are greatly reduced by 
using the HISQ action instead. 

The moderately heavy smearing used in the HISQ 
formalism complicates unquenched simulations. In Sec- 
tion VI, we discuss how these complications can be dealt 
with. Finally, in Section VII, we summarize our results 
and discuss their relevance to future work. Here we also 
present our first D physics results from the HISQ action. 

Much of our discussion is framed in terms of the 
"naive" discretization for quarks, which is conceptually 
equivalent to the staggered-quark discretization (see Ap- 
pendix [BJ . We usually find naive quarks to be more 
intuitive in analytic work than staggered quarks, while 
the latter formalism is definitely the more useful one for 
numerical work. We outline the formal connection be- 
tween the naive-quark and staggered-quark formalisms 
in a series of detailed Appendices. We use these results 
in our one-loop analysis of taste-exchange interactions, 
which is described in detail in Appendix IF1 

We used several sets of gluon configurations in this 
paper. The parameters for these various sets are sum- 
marized in Table [I] 
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a (fm) 


Gluon Action 


am s 


am u/d 


L/a 


UKQCD 


1/10 


Wilson 






16 


MILC 


1/8 


Liischer-Weisz 


0.05 


0.01 


20 


MILC 


1/11 


Liischer-Weisz 


0.03 


0.006 


28 



TABLE I: Gluon configurations used in this paper with in- 
formation about the collaboration that produced them, the 
gluon action used (the unimproved quenched Wilson action or 
the Symanzik-improved Liischer-Weisz action for full n/ = 3 
QCD), the s and u — d quark masses used in vacuum po- 
larization, and the spatial size of the lattice. The first set of 
configurations is described in |l3J, while the others are dis- 
cussed in detail in [Hll5|. 

II. SYMANZIK IMPROVEMENT FOR 
NAIVE/STAGGERED QUARKS 

A. Naive Quarks, Doubling, and Taste-Changing 
Interactions 

We begin our review of staggered quarks by examin- 
ing the formally equivalent "naive" discretization of the 
quark action (see Appendix |B|) : 

S = ^( x ) (7 • A(C7) + m ) i/>(x). (1) 
where A M is a discrete version of the covariant derivative, 

A M (£/)V(*) = 

— (U^(x) ip(x + ap) - Ufa - A) i>(x - ap)) , 

(2) 

Up(x) is the gluon link-field, a the lattice spacing, and 
too is the bare quark mass. The gamma matrices are 
hermitian, 

ll = 7m 7^ = 1 {7m>7<J = 2V> (3) 

where indices fj, and v run over 0. . . 3. A complete set of 
sixteen spinor matrices can be labeled by four-component 
vectors rip consisting of 0s and Is (i.e., rip £ 
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7n = n(7/,) n "- (4) 

A set of useful gamma-matrix properties is presented in 
Appendix [X] 

The naive-quark action has an exact "doubling" sym- 
metry under the transformation: 

rl>(x) -> ${x) = j^ p (-l) x "/ a i>(x) 

= 757 P exp(ix p ir/a)ip(x). (5) 

Thus any low energy-momentum mode, tp(x), of the the- 
ory is equivalent to another mode, ip{x), that has mo- 
mentum p p « ir/a, the maximum allowed on the lat- 
tice. This new mode is one of the "doublers" of the naive 



quark action. The doubling transformation can be ap- 
plied successively in two or more directions; the general 
transformation is 

ip(x) -> B c (x)ip(x) Jj(x) -» ${x)B\(x) (6) 

where 

B c (x) = 7? (-l) C ' x/0 

oc n (757 P ) Cp exp(ix ■ (ir/a), (7) 
p 

and £ is a vector with £ M £ Z2, while £ is "conjugate" 
to £ (see Appendix 1X1) : 

C = E Cu mod 2 (8) 

Consequently there are 15 doublers in all (in four di- 
mensions), which we label with the fifteen different 
nonzero £'s. 

As a consequence of the doubling symmetry, the stan- 
dard low-energy mode and the fifteen doubler modes 
must be interpreted as sixteen equivalent flavors or 
"tastes" of quark. (The sixteen tastes are reduced to four 
by staggering the quark field; see Appendix [Bj) This 
unusual implementation of quark tastes has surprising 
consequences. Most striking is that a low-energy quark 
that absorbs momentum close to C 71 "/ : f° r one of the 
fifteen £'s, is not driven far off energy-shell. Rather it is 
turned into a low-energy quark of another taste. Thus 
the simplest process by which a quark changes taste is 
the emission of a single gluon with momentum q sa £ir/a. 
This gluon is highly virtual, and therefore it must imme- 
diately be reabsorbed by another quark, whose taste will 
also change (see Fig.^l. Taste changes necessarily involve 
highly virtual gluons, and so are both suppressed (by 
a 2 ) and perturbative for typical lattice spacings. Onc- 
gluon exchange, with gluon momentum q « C 71 "/ ; is 
the dominant flavor-changing interaction since it is low- 
est order in a s (£7r/a) and involves only 4 external quark 
lines. (Processes with more quark lines are suppressed by 
additional powers of (ap) 3 where p is a typical external 
momentum.) This observation is crucial when trying to 
improve naive quarks by removing finite-a errors, as we 
discuss below (see also p.l9llTo|). 

Sixteen tastes of quark from a single quark field is fif- 
teen too many. In the absence of taste exchange, factors 
of 1/16 are easily inserted into simulations to remove the 
extra copies. In particular, quark vacuum polarization is 
corrected by replacing the quark determinant in the path 
integral by its 1/16 root, 

det(7 ■ A + to ) -> det(7 • A + m ) 1/16 , (9) 

or, equivalently, by multiplying the contribution from 
each quark loop by 1/16. Recent empirical studies of 
the spectrum of 7 • A, and improved versions of it, show 
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(ir/a 



sequently, while taste exchange can lead to such anoma- 
lies as unitarity violations, these are suppressed in the 
continuum limit. They are further, and more efficiently, 
suppressed through Symanzik improvement of the action, 
as we now discuss. 




FIG. 1: The leading tree-level taste-exchange interaction, 
which involves the exchange of a gluon with momentum Qir/a 
where each f M is or 1 but £ 2 7^ 0. 



how eigenvalues cluster into increasingly degenerate mul- 
tiplets of sixteen as the lattice spacing vanishes, which is 
necessary if this procedure is to be useful for quark vac- 
uum polarization [T^ . The treatment of valence quarks 
is illustrated in Appendix IU1 

It is easy to see how the 1/16 root corrects for taste 
in particular situations. For example, consider the vac- 
uum polarization of a single pion in a theory of only one 
massless quark flavor. Ignoring taste exchange, the dom- 
inant infrared contributions come from the quark dia- 
grams shown in Fig. [21 (gluons, to all orders, are im- 
plicit in these diagrams). The first diagram, with its 
double quark loop, has contributions from 16 2 identi- 
cal massless pions, corresponding to configurations where 
the quark and antiquark each carries a momentum close 
to an integer multiple of ir/a (and is therefore close to be- 
ing on shell). This diagram is multiplied by 1/16 2 since 
there are two quark loops, thereby giving the contribu- 
tion of a single massless pion. The second diagram in- 
volves quark annihilation into gluons, and has contribu- 
tions from 16 identical massless pions, corresponding to 
configurations where the quark and antiquark each car- 
ries a momentum close to the same C 71 "/ ( so they can 
annihilate into low- momentum gluons). This diagram is 
multiplied by 1/16 since there is only one quark loop, 
thereby giving the contribution again of a single mass- 
less pion but now with an insertion from the gluon decay. 
Inserting additional annihilation kernels results in a ge- 
ometric series of insertions in the propagator of a single 
massless pion that shifts the pion's mass away from zero, 
as expected in a U(l) flavor theory. 

Such patterns are disturbed by taste-exchange inter- 
actions like that in Fig. ^ For example, these interac- 
tions cause mixing between the different tastes of pion 
in Fig. |3 This mixing lifts the degeneracy in the pion 
masses so that different tastes of pion are only approx- 
imately equivalent. Consequently the "1/16-root rule" 
gives results that correspond only approximately to a 
single pion. Taste-changing interactions are suppressed 
by a 2 , however, and detailed analyses, using chiral per- 
turbation theory, demonstrate that the impact on physi- 
cal quantities is therefore also suppressed by a 2 |l8j . Con- 



B. Tree-Level Symanzik Improvement 

The discretization errors in the naive-quark action 
come from two sources. The more conventional of these 
corrects the finite-difference approximation to the deriva- 
tives in the action: one replaces :6j 



A, 



" A ^ 

— A J 

6 " 



(10) 



in the naive-quark action (Eq. JIJ). The a 2 correction is 
often referred to as the "Naik term." 

Less conventional is a correction to remove leading or- 
der taste-exchange interactions As discussed 
in the previous section, these interactions result from the 
exchange of single gluons carrying momenta close to Qit /a 
for one of the fifteen non-zero (s (^ G Z2). Since these 
gluons are highly virtual, such interactions are effectively 
the same as four-quark contact interactions and could be 
canceled by adding four-quark operators to the quark 
action. These operators affect physical results in 0(a 2 ) 
since they have dimension six. A simpler alternative to 
four-quark operators is to modify the gluon-quark ver- 
tex, ipj^U^ip + • • • , in the original action by introducing 
a form factor fp(q) that vanishes for (taste-changing) 
gluons with momenta q = C 7T / a f° r each of the fifteen 
nonzero £s. In fact the form factor for direction p, need 
not vanish when £ M = 1 since the original interaction 
already vanishes in that case. Consequently we want a 
form factor where 



1 for q — > 

for q^(Tr/a where C 2 ^ 0, C M = 0. 



(11) 




We can introduce such a form factor by replacing the 
link operator Up{x) in the action with J-pUp(x) where 
smearing operator Tp is defined by 



(12) 



(2) 

and 5 P approximates a covariant second derivative when 
acting on link fields: 

4 2) Up{x) = i (Up(x) Up(x + ap) U\{x + ap) 
- 2Up(x) 

+ Wpix — ap) Up(x — ap) U p (x — ap + a/t)) . 

(13) 
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FIG. 2: Infrared contributions from (naive lattice) quark vacuum polarization corresponding to a pion vacuum polarization 
loop in a simulation of one flavor. The sums are over all infrared sectors, where quarks have momenta near p M = (^ir/a for one 
of the sixteen £s consisting of just Os and Is. The factors of 1/16 are from the 1/16 rule (Eq. ©) and cancel the sums. Gluons, 
to all orders, are implicit in these diagrams. 



This works because 8 p 2 ^ « —A/a 2 (and J-^ vanishes) when 
acting on a link field that carries momentum q p « n/a. 
This kind of link smearing is referred to as "Fat 7" smear- 
ing in pfl. 

Smearing the links with removes the leading 
0(a 2 ) taste-exchange interactions, but introduces new 
Q(a ) e rrors. These can be removed by replacing J-^ 
with [13 



jrASQTAD _ j: _ \^ fl2 (^p) 2 



(14) 



where 5 p approximates a covariant first derivative: 

5 P U M (x) = - (U p (x) U^(x + ap) U%(x + ap) 

— UUx — ap) Up{x — ap) U p (x — ap + ap,)) 



(15) 



The new term has no effect on taste exchange but (ob- 
viously) cancels the 0(a 2 ) part of T^. Correcting the 
derivative, as in Eq. (|l(Jfl . and replacing links by a 2 - 
accurate smeared links removes all tree-level O (a 2 ) errors 
in the naive-quark action. The result is the widely used 
"ASQTAD" action [13, 



where, in the first difference operator, 
V,(x)^Tf^U,(a 




In practice, operator Vn is usually tadpole improved |16| : 
in fact, however, tadpole improvement is not needed 
when links are smeared and reunitarized 19, 201. 



C. c Quarks 

The tree-level discretization errors in the ASQTAD ac- 
tion are 0((ap^) 4 ), which are negligible (< 1%) for light 
quarks. When applied to c quarks, however, these errors 
will be larger. The most important errors in this case 
are associated with the quark's energy, rather than its 
3-momentum, since c quarks are typically nonrelativistic 
in ips, Ds and other systems of current interest. Conse- 
quently E « m p and the largest errors are 0((am) 4 ) 
or 6% at current lattice spacings (am c rs 0.5). Such er- 
rors show up, for example, in the tree-level dispersion 
relation of the quark where: 



c 2 (0) ee lim 



E 2 (p) 



= l + l( am )4 + I (flTO) 6 + 



(18) 



This dominant (tree-level) error can be removed by 
retuning the coefficient of the a 2 7^A^ (Naik) term in the 
action: 

]T ?(*) (£7, (X - £(l + e)A») +m j V(*), 



where parameter e has expansion 



27 . x o , 327 . .a 

— (am) H (am) 

40 v ; 1120 ^ ; 

153607 



5843 
53760 



{amf 



3942400 



{amf 



(19) 



(20) 



With this choice, c 2 = l + 0((am) 12 ). Only the first term 
(17) in e's expansion is essential for removing (am) errors; 



the remaining terms remove errors to higher orders in am 
but leading order in v/c, the quark's typical velocity. In 
practice, terms beyond the first one or two are negligible 
though trivial to include. 
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Note that the bare mass Too in the quark action is 
related to the tree-level quark mass by 



III. SYMANZIK IMPROVEMENT AT 
ONE-LOOP ORDER 



m n — ml 1 H (am) (am) + 

1 80 V ' 2240 1 ' 



(21) 



for this choice of e. This formula is useful for extracting 
the corrected bare quark mass, to, from a tuned value 
of Too- 

To examine the tree-level errors in this modified ASQ- 
TAD action, it is useful to make a nonrelativistic expan- 
sion since, as we indicated, heavy quarks are generally 
nonrelativistic in the systems of interest. This expansion 
is easy because the action is identical in form to the con- 
tinuum Dirac equation, and consequently has the same 
nonrelativistic expansion but with 



One-loop errors, which are 0(a s (ap) 2 ), should be no 
more than a percent or so for light quarks and current 
lattice spacings. Taste-changing errors, however, have 
generally been much larger than expected. These errors 
also bear directly on the validity of the 1/16-root trick 
for accounting for taste in the vacuum polarization. For 
this reason it is worth trying to further suppress taste- 
exchange errors beyond what has been accomplished at 
tree level. We show how to do this in this section. 

The only other one-loop effects that are important at 
the level of a few percent are a s (am) 2 for heavy quarks, 
and, in particular, the charm quark. We show how to 
correct for these errors in the Section lill Bl 
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A. Light Quarks 



At tree level, 

— > smh(aD t — am) /a, 

when we replace 

ip — ► exp(— mt) tp, 



(23) 



(24) 



where to is the quark mass. This redefinition removes the 
rest mass from all energies. Expanding in powers of 1/to, 
we then obtain the tree-level nonrelativistic expansion of 
the modified ASQTAD action: 



NR 



D, 



D 2 

2m 



D 4 



ck 



ga ■ B 



ig[D-,E] 

+ce — o — 5 CE 

8to 2 



8to 3 2m 
ga ■ (D x E - E x D) 

8m 2 



where 



ck 
ce 



— (am) H (am) — 

10 y ' 280 v ' 

A ( am f + JL (am) 6 
80 V ; 2240 V ; 



4>NR 

(25) 

(26) 
(27) 



The coefficients of all terms through order 1 /to have no 
am errors at tree level. In particular the a ■ B term is 
correct at tree level once to is tuned (nonperturbatively) . 
The coefficients of the remaining terms shown are off by 
less than 6% and 1% for charm quarks on current lat- 
tices (am ps 1/2), which should lead to errors of order 
(am) 4 (v/c) 2 w 2% or less in charmonium analyses and 
less than 0.5% for D physics. (Errors should actually be 
considerably smaller than this because of dimensionless 
factors like the 1/8 in the D 4 term.) Other tree-level er- 
rors, which are order (ap) 4 and e(ap) 2 /6, should also be 
less than 1% for a w 0.1 fm. 



One-loop taste exchange comes only from the diagrams 
shown in Fig. |3| other one-loop diagrams vanish because 
of the tree-level improvements. Again the gluons in these 
diagrams transfer momenta of order Cit/a, and so are 
highly virtual. Consequently, the correction terms in the 
corrected action that cancel these will involve current- 
current interactions that conserve overall taste (because 
of momentum conservation). There are 28 such terms 
in the massless limit for the staggered- quark action (see 
Appendix : 
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J, 
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(5 M ) 
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(m) 



(5) 



+ (color octet Js with d$ -> , 
where sums over indices are implicit and v ^ /i. The cur- 

(s) 

rents Jn in the first 14 operators are color-singlets with 
taste s and spinor structure n, while the last 14 operators 
are the same but with color octet currents. The precise 
definitions of these currents are given in Appendix lEl 

(s) ~(s) 

We have computed the coupling constants d n and d„ 
for the tree-level improved theory; the results are pre- 
sented in Table [n] (see also Appendix • Armed with 
these one could incorporate A£ contac t into simulations, 
removing all one-loop taste exchange. This procedure is 
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FIG. 3: The only one-loop diagrams that result in qq — > qq taste exchange for ASQTAD quarks, 
external quark lines are suppressed by additional powers of a 2 . 



Diagrams with additional 



complicated to implement, however. A simpler procedure 
is suggested by the results in That analysis shows 
that repeated smearing of the links further reduces the 
mass splittings between pions of different taste. 

This result is not surprising given the perturbative ori- 
gins of taste exchange. The dominant taste-exchange 
interactions (in ASQTAD) come from the one-loop di- 
agrams in Fig. |3 The largest contributions to these di- 
agrams come from large loop momenta, k — 0(Cit /2a). 
Smearing the links introduces form factors that suppress 
high momenta, leaving low momenta unchanged, and 
therefore should suppress this sort of one-loop correction. 
The links in the ASQTAD action are smeared, but the op- 
erator 

^ASQTAD which 

smears links in the y, direction, 
introduces additional links in orthogonal directions (to 
preserve gauge invariance) that are not smeared. Smear- 
ing the links multiple times guarantees that all links are 
smeared, and high momenta are suppressed. 

There are three problems associated with multiple 
smearing. One is that it introduces new 0(a 2 ) dis- 
cretization errors, and these errors grow with additional 
smearing. The action presented in suffers from 
this problem. The a 2 errors can be avoided by using 
an a 2 -improved smearing operator, such as ^ AS Q TAD 
(Eq. IjMfl ). Replacing smearing operator 



jrASQTAD 



^-ASQTAD ^-ASQTAD 



(29) 



in the quark action, for example, smears all links but 
does not introduce new tree-level a 2 errors. 

The second problem with multiple smearings is that 
they replace a single link in the naive action by a sum of 
a very large number of products of links. This explosion 
in the number of terms does not affect single-gluon ver- 
tices on quark lines, by design, but leads to a v^V-growth 
in the size of two-gluon vertices where N is the number 



of terms in the sum. The v iV-growth enhances one- loop 
diagrams with two-gluon vertices, canceling much of the 
benefit obtained from smearing. This problem is reme- 
died by reunitarizing the smeared link operator: that is 
we replace 



^-ASQTAD 



^•ASQTAD jj ^ASQTAD 



(30) 



where operator hi unitarizes whatever it acts on. A 
smeared link that has been unitarized is bounded (by 
unity) and so cannot suffer from the VN problem. Al- 
though we verified that it is unnecessary to make the 
unitarized link into an SU(3) matrix — simple unitariza- 
tion is adequate — the results given in this paper use links 
that are projected back onto ST/(3). Reunitarizing has 
no effect on the single-gluon vertex, and therefore does 
not introduce new 0(a 2 ) errors. 

The doubly smeared operator is simplified if we rear- 
range it as follows 



a 2 {5 p f 



UT, 



(31) 



where the entire correction for a 2 errors is moved to the 
outermost smearing. Our new "HISQ" discretization of 
the quark action is therefore 



J2 1>(x) ( 7 -2? HISQ +m) 

X 

where 

pHISQ = A ^ (w) _^ (1 + e)A 3 (x) 

and now, in the first difference operator, 



(32) 

(33) 
(34) 
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Unimproved Gluons Improved Gluons 
ASQTAD HISQ HYP ASQTAD 



Octet: 





1.41 


0.19 


0.02 


0.77 


<c 


0.46 


0.00 


0.08 


0.25 




0.23 


0.01 


0.00 


0.16 


"-"[IIS 


0.38 


0.02 


0.01 


0.27 


<gr> 


0.34 


0.08 


0.04 


0.19 




0.35 


0.03 


0.02 


0.19 




0.08 


0.00 


0.01 


0.05 




0.20 


0.00 


0.01 


0.13 


d M 


0.20 


0.01 


0.02 


0.11 




0.31 


0.01 


0.01 


0.17 




0.06 


0.00 


0.00 


0.04 


d M 
"5 


0.07 


0.00 


0.01 


0.04 


d (5) 
u 5m 


0.17 


0.00 


0.01 


0.09 


4 5) 


0.06 


0.00 


0.00 


0.04 



Singlet: 










j(5/u) 
"5 


0.76 


0.10 


0.01 


0.41 


d^ 


0.12 


0.00 


0.00 


0.09 


d { r ] 


0.18 


0.05 


0.02 


0.10 




0.04 


0.00 


0.00 


0.03 


d (v) 


0.11 


0.01 


0.01 


0.06 




0.03 


0.00 


0.00 


0.02 


d (5) 


0.09 


0.00 


0.00 


0.05 



Avg d,d 0.23 0.02 0.02 0.13 
a 2 Aml 0.023(1) 0.011(1) 0.009(1) 



TABLE II: 0(a s ) coefficients for the couplings in £ C ontact 
(Eq. 1281 ) for three different light-quark actions and massless 
quarks. Terms not listed have zero coefficient in this order. 
Results are given for the unimproved (Wilson) gluon action, 
and, in the ASQTAD case, also for the standard Liischer- 
Weisz action. Results for the average coefficient are also given, 
as are the splittings Amjj. between the 3-link pion and the 
Goldstone pion. The pion masses are from simulations with 
uncorrected gluons at 1/a » 1/10 fm with no vacuum polar- 
ization. 



To examine this issue, we recomputed the one-loop co- 
efficients for all of the leading taste-exchange correction 
terms in the action (Eq. I|28l) ) using the HISQ action, 
as well as the HYP action from The results are 

derived in Appendix [0 and presented in Table ^ The 
HISQ coefficients are on average an order of magnitude 
smaller than the ASQTAD coefficients. Consequently 
the residual taste-exchange one-loop corrections in the 
HISQ action are probably smaller than two-loop correc- 
tions and possibly also dimension-8 taste-exchange oper- 
ators at current lattice spacings. The HYP action gives 
coefficients similar in size to HISQ. 



B. c Quarks 

Many new taste-preserving operators would need to 
be added to the HISQ action in order to remove all or- 
der a s a 2 errors. The leading operators are listed in Ta- 
ble [ffi] None of these is relevant at the few percent level 
for light quarks, but terms that enter in relative order 
a s {am c ) 2 might change answers by as much as 5-10% 
for c quarks on current lattices, where am c ~ 0.5 and 
a s « 1/3. High-precision work requires that these errors 
be removed. 

Nonrelativistic expansions of the various operators in 
Table IIHI show that only the first, the Naik term, can 
cause errors in order a s (am c ) 2 . All of the others result 
in harmless renormalizations or are suppressed by addi- 
tional powers of v/c, and so contribute only at the level 
of 2-3% or less for and less than 1% for Ds. We can 
remove all a s (am c ) 2 errors by including radiative correc- 
tions in the e parameter that multiplies the Naik term 
(Eq. Q2J}). 

There are two ways to compute the radiative correction 
to e. One is nonperturbative: e is adjusted until the 
relativistic dispersion relation, 

x E 2 (p)-m 2 , x 

c 2 (p) = —^L = 1, (36) 

P 



while in the second, 

X^^UT^U^x). (35) 

The third problem with smearing is that previous work 
has only demonstrated its impact on the mass splittings 
between pions with different tastes. For small masses, 
there are only five degenerate multiplets of pion. Con- 
sequently showing that smearing reduces the splittings 
between the pion multiplets does not guarantee that all 
28 taste-exchange terms in the quark action (Eq. (J2HJ) 
are suppressed; and, therefore, it does not guarantee that 
taste-exchange effects are generally suppressed (other 
than in the pion mass splittings). 



computed in a meson simulation is valid for all low three- 
momenta p. 

A simpler procedure is to compute the one-loop cor- 
rection to e using perturbation theory — by requiring, for 
example, the correct dispersion relation for a quark in 
one-loop order. The result would have the form 

27 

e = e ± a s - —(am) 2 + 0(a 2 ,(am) 4 ) (37) 

where ej depends upon (am) 2 . As we will show, ei turns 
out to be negligibly small for the HISQ action (but not 
for ASQTAD). 
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Operator 


ip Physics 


D Physics 


a 2 ^D^7 M V 


a s (am) 2 


a s (am) 2 


a 2 mipa ■ gFtp 


a s (v/c) 2 (am) 2 


a 3 (v/c)(am) 2 


a 2: ipD 2 D ■ -yip 






a 2 %p{D ■ y) 3 ip 






a mipD tp 






a 2 ipa ■ gF D ■ jtp 






a 2 ipD ■ gF ■ ftp 


a s (v/c) 2 (am) 2 


a s (v/c) 2 (am) 2 



a 2 (ipyip) 2 
a 2 (V'775V') 2 



TABLE III: Relative errors associated with the leading (for 
charm quarks) taste-preserving operators that enter the HISQ 
action in one-loop order. Note that a s ~ 1/3, (am) 2 ~ 1/4, 
and (v/c) 2 ~ 1/3 for ips (1/10 for Ds) for lattice spacings 
of order 0.1 fm. The errors listed are relative to the binding 
energy which is of order 500 MeV for both tps and Ds. Some 
of the operators are related to others in the same grouping by 
the equations of motion, and so are redundant. 



IV. APPLICATIONS: BOUNDING TASTE 
EXCHANGE 

We tested our perturbative analysis of the suppression 
of taste exchange by computing the pion mass taste split- 
tings for ASQTAD, HISQ, and a few other variations |21| . 
Some of these results are shown in Table[H]where we list 
the mass splitting AmJ between the 3-link pions and 
the 0-link (Goldstone) pion from a quenched simulation 
with unimproved (Wilson) glue at a = 1/10 fm and quark 
mass am = 0.03 [13J. As expected, small mass splittings 
are correlated with small values for the coefficients of the 
one-loop taste-exchange interactions (Eq. i|28|) ). These 
results confirm: a) that taste-exchange interactions are 
perturbative in origin; b) that these interactions will be 
suppressed by a factor of order a s ~ 1/3 in the HISQ 
action relative to the ASQTAD action at current lattice 
spacings; and c) the HISQ action will be more accurate 
than other current competitors because it has negligible 
one-loop taste exchange and no tree- level 0(a 2 ) errors. 

We have also compared meson splittings for ASQ- 
TAD and HISQ valence quarks in full QCD simula- 
tions with rif = 3 light-quark (ASQTAD) vacuum po- 
larization. We find that the pseudoscalar splittings are 
3.6(5) times smaller with HISQ than with ASQTAD on 
lattices with a — 1/8 fm, when both valence quarks are 
s quarks. The splittings are 3.1(7) times smaller on our 
1/11 fm lattices. 

Any taste-exchange effect will be 3-4 times smaller in 
HISQ, and so measuring the ASQTAD-HISQ difference 
in any quantity allows us to bound taste-exchange errors 
in each formalism. To illustrate how such comparisons 
are done, we have computed the mass of the <p meson 
and the pseudovector decay constant f rh of the rj s meson 
using ASQTAD and HISQ valence quarks for each of our 
two lattice spacings. 



The <p is unstable and so its mass in a simulation is sen- 
sitive to the volume of the lattice (it is not "gold-plated" 
in the sense of [lj). We have not corrected for this vol- 
ume dependence, but have simulated using a fixed phys- 
ical volume. Consequently results for the mass should 
extrapolate to the same value as a — > 0, but this value 
may be a few tens of MeV larger than the physical mass. 
The rj s is an easy-to-analyze substitute for a 7r; it is a 
h ss meson, but without valence-quark annihilation, 
and therefore hypothetical. We tuned the s mass so that 
the r) s masses in each case were all identical (696 MeV) 
to within ±lMeV (statistical/fitting errors only). We 
evaluated the rj s decay constant by computing the ma- 
trix element of the pseudoscalar density and multiplying 
by the bare quark mass (see Q , for example) . 

Our results, shown in Table IIVI suggest significant 
taste-exchange errors for ASQTAD, particularly in the 
1/8 fm simulation. Looking first at the decay con- 
stant f r)s , the ASQTAD result on the coarse lattice 
is 6(3)% larger than the ASQTAD result on the fine lat- 
tice, where taste-exchange effects should be at least a 
factor of two smaller. The uncertainty in the discrep- 
ancy is large here because of uncertainties in the val- 
ues of the lattice spacings. Far more compelling is the 
4.6(6)% discrepancy between ASQTAD and HISQ on the 
coarse lattice, which can be measured far more accurately 
than differences between simulations with different lat- 
tice spacings. This discrepancy is reduced to 1.9(4)% 
on the fine lattice, as expected since a 2 is half as large. 
Since errors in HISQ should be 3-4 times smaller than in 
ASQTAD, these discrepancies imply that taste-exchange 
errors in the decay constant could be as large as 5-6% 
for ASQTAD but only 1-2% for HISQ when a = 1/8 fm. 
The errors are half that size when a = 1/11 fm. Results 
from both formalisms are consistent with an extrapolated 
value of 180 (4) MeV. 

Results for the <f> mass are very similar. These give a fi- 
nal mass of 1.05(2) GeV, which is 30(20) MeV above the 
mass 1.019 GeV from experiment (and consistent with 
expectations for finite- volume effects). These analyses 
taken together indicate that HISQ is delivering 1-2% pre- 
cision already on the coarse lattice, while a significantly 
smaller lattice spacing (and much more costly simulation) 
is needed to achieve similar precision from ASQTAD. 
HISQ a 2 errors should be less than 1% on the 1/11 fm 
lattice for light quarks. 

Our tests are only partial because ASQTAD quarks 
were used for the vacuum polarization in both the HISQ 
and ASQTAD analyses, but the effects of changes in 
the vacuum polarization are typically 3-5 times smaller 
than the effects caused by the same changes in the va- 
lence quarks. The agreement between our HISQ results 
from different lattice spacings indicates that finite-a er- 
rors from the vacuum polarization are not large for these 
lattice spacings. 

In past work, taste-exchange errors have been esti- 
mated by comparing results from different lattice spac- 
ings, relying upon the fact that these errors vanish as 
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a (fm) ASQTAD HISQ ASQTAD/HISQ 



f v , (MeV) 1/8 196(4) 187(4) 1.046(6) ASQTAD (1/11 fm): 

1/11 185(3) 182(3) 1.019 (4) 

m (GeV) 1/8 1.119(24) 1.076(23) 1.040(9) 

1/11 1.057(18) 1.052(16) 1.005(6) HISQ (1/11 fm): 



TABLE IV: Simulation results for decay constant f V3 and 
<f> mass using either ASQTAD or HISQ for the valence quarks 
from full QCD simulations (with n/ = 3 ASQTAD quarks). 
Lattice spacing errors cancel in the ratio ASQTAD to HISQ 
computed on the same gluon configurations; thus errors are 

much smaller for the ratio. HISQ (1/8 fm): 



a — > 0. This example illustrates how the same errors can 
be reliably estimated by comparing ASQTAD with HISQ 
results from simulations using only a single lattice spac- 
ing. This approach to estimating taste-exchange errors 
is very efficient. It is almost certainly the simplest way 
to quantify uncertainties about taste exchange and the 
validity of the 1/16-root trick for vacuum polarization. 



c 


am 


(p/Pmin) 2 


c 2 (p) 





0.38 


2 


0.962 ( 9) 


0.3 


0.38 


2 


1.020 (11) 


0.4 


0.38 


2 


1.036 (11) 





0.43 


2 


1.029 (11) 


-0.115 


0.43 


1 


0.985 (16) 


-0.115 


0.43 


2 


0.992 (13) 


-0.115 


0.43 


3 


1.014(15) 


-0.115 


0.43 


4 


0.991 (11) 





0.67 


1 


1.190 (20) 


-1 


0.67 


1 


0.560 (10) 


-0.35 


0.67 


1 


0.904 (15) 


-0.28 


0.67 


1 


0.950(15) 


-0.21 


0.66 


1 


1.008 (13) 


-0.21 


0.66 


2 


1.017(10) 


-0.21 


0.66 


3 


1.019(12) 


-0.21 


0.66 


4 


1.007 ( 7) 



V. APPLICATIONS: c QUARKS AND 
CHARMONIUM 

As argued in Sections III CI and IIII Bl we expect the 
ASQTAD action and, especially, the HISQ action to work 
well for c quarks even though am c is typically of order 
0.4 or larger on current lattices. To achieve high pre- 
cision (few percent or better), we must tune the Naik 
term's renormalization parameter e (Eq. I|19(l ). Here we 
do this nonperturbatively by computing the speed of light 
squared, c 2 (p) (Eq. (|18|) ). for the r\ c in simulations with 
various values of e and tuning e until c 2 = 1. Our results 
are summarized in Table Ivl 

The ASQTAD results on the fine lattices show only 
small errors in c 2 even with e = 0. Values for different 
es are plotted in Fig. together with an interpolating 
curve. These data indicate that the optimal choice is e — 
0.19(5) for this lattice spacing and mass. The tree-level 
prediction for e, from Eq. (|20[1 . is —0.10, which indicates 
that radiative corrections in e are of order 1 x a s ~ 0.3, 
as expected. 

The HISQ results show even smaller c 2 errors on the 
fine lattices. Tuning to e = —0.1155, the tree-level value 
given by Eq. (|20[l for am c — 0.43, removes all errors 
in c 2 at the level of 1%. This suggests that one-loop 
and higher-order radiative corrections in e are negligi- 
ble for HISQ compared with the tree-level corrections. 
The dominance of tree-level contributions is confirmed 
by our HISQ analysis using the coarser lattice spacing, 
where am c w 0.66. Here, we found the optimal value 
e = —0.22(3) again by tuning e until c 2 = 1 for a low- 
momentum j] c . (We overshot slightly and did simulations 
for e = —0.21 rather than —0.22.) Our tuned e compares 
quite well with the tree-level prediction of —0.246 from 



TABLE V: The speed of light squared, c 2 (p), computed 
from the Tfc's dispersion relation in simulations using differ- 
ent quark actions. Results are given for different values of 
the Naik term's renormalization parameter e (Eq. 1191 '). the 
lattice spacings a, and the meson's three-momenta p (in units 
of the smallest momentum on the lattice which was roughly 
0.5 GeV in each case). 



Eq. i|20|) . Consequently it is quite likely that the tree- 
level formula is sufficiently accurate for most practical 
applications today. 

Note that setting e = — 1 cancels out the Naik term 
completely. Tree-level errors are then order a 2 rather 
than order a 4 , as in the HISQ action. Table IVI shows 
that these a 2 errors cause c 2 to be off by almost a factor 
of 2. This example underscores the importance of using 
a 2 -improved actions in high-precision work. 

On the coarse lattice, we use the c 2 (p) with p = p m i n , 
the smallest nonzero momentum on our lattice, to tune e. 
It is important to verify that the tuned action gives the 
correct dispersion relation for other momenta as well. 
The data in Fig. demonstrate that errors are less than 
a couple percent for meson momenta out to 1 GeV even 
on the coarse lattice where am c = 0.66. Indeed, the er- 
rors would probably have been smaller (< 1%) had we 
tuned e a little more accurately. 

We examined the spectrum of the ip meson family to 
further test the precision of our formalism. We used the 
1/11 fm lattices where am c w 0.4, and therefore we ex- 
pect errors of order 1-2% of the binding energy, or 5- 
10 MeV (see Section fill B|) . A particularly sensitive test 
is the hyperfine mass splitting between the ip and rj c . We 
tuned the bare c mass until the mass of the local (and 
lightest) t] c in our simulation agreed with experiment. 
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0.95 


- cj) e = 0.19(5) 

1 , 1 







0.2 



0.4 



FIG. 4: The speed of light squared, c 2 (p), computed for an 
77c meson using the ASQTAD action with p = (p m in,Pmm, 0) 
and different Naik-term renormalizations e, where p m in is the 
smallest lattice momentum possible (approximately 500 MeV 
here). The correct value, e = 0.19(5), occurs at the intersec- 
tion of the interpolating line with the line c 2 = 1, as shown. 
Here am c = .38. 



other than pseudoscalars are much less sensitive to taste- 
exchange effects. 

The ip — rjc hyperfine splitting, for 0-link mesons in 
each case, is 109(3) MeV in our 1/11 fm simulation with 
tuned e (= —0.115). The 3% uncertainty is almost en- 
tirely due to tuning uncertainties in the lattice spacings 
since these uncertainties enter twice: once for convert- 
ing the splitting from lattice to physical units, and once 
through uncertainties in the c mass, which are themselves 
controlled by uncertainties in a -1 |12|. Our HISQ result 
is somewhat smaller than the current experimental result 
of 117(1) MeV 23j, but it needs three further corrections. 
(It is also worth noting that the current Particle Data 
Group average of 117 MeV is somewhat larger than the 
most recent experiments which find values in the range 
113-115 MeV with uncertainties of 1-2 MeV; see pl|.) 

The first correction comes from the operators in Ta- 
ble 11111 that must be added to the HISQ lagrangian in 
order to remove further discretization errors. Of these 
the most important for the hyperfine splitting is 



1.1 



0.9 




SCms = Chfs a 2 m c ipa ■ Ftp 



(38) 



0.25 



0.5 0.75 
> 2 (GeV 2 ) 



FIG. 5: The speed of light squared, c (p ), for rj c mesons 
at different momenta on a lattice where am c ~ 0.66 for the 
HISQ action with e — —0.21. A comparison with other results 
from Table IV! suggests that choosing e = —0.22 would move 
all points down onto the line c 2 = 1 to within errors. 



Since there are no other free parameters in our action, 
the simulation then predicts a mass for the ip. As is 
clear from Fig. the ip mass in the simulation is quite 
accurate — certainly within the 5-10 MeV we expected. 

The data shown in this last figure also bound taste- 
exchange errors. As shown there, different tastes of the rj c 
have different masses, because of taste-exchange interac- 
tions. The maximum spread in the masses for the HISQ 
action, however, is only 9 MeV. This is much smaller than 
the spread from the ASQTAD action, where the split- 
ting between, for example, the 0-link and 1-temporal- 
link pion masses is 40 MeV (compared to only 3 MeV for 
the same splitting with HISQ). (Uncorrected staggered 
quarks also show very large splittings Note also 

that the spread in the ips is three times smaller than the 
spread in the rj c s. This is typical; the masses of mesons 



where Chfs = 0(a s ). This will affect the hyperfine split- 
ting in relative order a s (am c ) 2 , or at the level of ±5 MeV. 
The coefficient Chf s is readily computed in perturbation 
theory and this calculation is underway. 

The second correction comes from residual taste- 
exchange interactions, which from the data in Fig. 
could be of order a few MeV. Both this error and that 
corrected by SC^h are approximately proportional to a 2 . 
So we can estimate them (together) by comparing to a 
calculation with a different lattice spacing. We repeated 
our hyperfine splitting analysis on the 1/8 fm lattice us- 
ing the tuned e (= —.21). The a 2 errors should be a little 
more than twice as large on the coarser lattice since a 2 
is almost exactly twice as large. We obtained a splitting 
of 110(3) MeV on the coarse lattice, which is essentially 
identical to the 109(3) MeV we obtained on the fine lat- 
tice. Combining these results, together with our a priori 
expectation of ±5 MeV errors from a 2 corrections, we ob- 
tain an a 2 corrected hyperfine splitting is 109(5) MeV. 

The third correction is due to the fact that our simula- 
tion does not include effects from the annihilation of the 
valence cc quarks into two or more gluons. Such anni- 
hilations are responsible for small shifts in the r\ c and ip 
masses, as well as for the (non-electromagnetic) hadronic 
decay rate of each meson. The dominant contribution 
comes from cc — > gg and affects only the r\ c . The shift in 
the ?7 C energy is proportional to the perturbative ampli- 
tude for cc — > gg — * cc at threshold [24J and therefore [2{| 



AE nc ~iT(r] c —>■ hadrons)/2 



ocln(2) - l-i- + 0(a s ). 



(39) 



This result implies that the leading correction to the r\ c 
mass due to cc annihilation can be computed from the 
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FIG. 6: Masses for different tastes of the r\ c and V using 
HISQ and ASQTAD c quarks with e = on a 1/11 fm lattice 
(am c »s 0.4). For HISQ, ?7 c s are given (from left to right) 
for the 0, 1, 2, and 3 spatial link tastes, without and then 
with a temporal link; only the spatial splittings are given 
for the ip. For ASQTAD, r) c results are given for only the 
O-link and 1-temporal-link operators, and i\> results for only 
the 0-link operator. The dashed lines indicate the results 
from experiment. Error bars are of order the size of the plot 
symbols. 



hadronic width: 



> 

CP 

O 



3.8 



3.6 



3.4 



3.2 




FIG. 7: Masses for different excitations of the ip meson from 
a simulation at lattice spacing 1/11 fm. The dashed lines in- 
dicate the results from experiment. 



AE Vc = T(r; c -> hadrons) f — - + 0{a s 
= -2.4(8) MeV 



(40) 



where we use experimental results from |23j . This correc- 
tion increases our theoretical value for the ip-r] c splitting 
to 111(5) MeV, which agrees well with experiment. The 
sea-quark masses in our simulations are not quite cor- 
rect, but both theoretical expectations and experience 
with previous simulations indicate that this has negligi- 
ble effect on this hyperfine splitting. 

This level of precision would be impossible using our 
ASQTAD results because the taste-exchange errors are 
tens of MeV — much larger than the ±5 MeV expected 
from other a 2 errors. This example confirms that taste- 
exchange errors are likely the dominant source of a 2 er- 
rors in the ASQTAD formalism. With HISQ, on the 
other hand, taste-exchange errors have been suppressed 
to a level commensurate with other errors. 

Finally we also computed the masses of some radially 
and orbitally excited states in the ip family using HISQ, 
although not as accurately as the ground-state masses. 
High-precision determinations of excited-state masses re- 
quire careful design of the meson sources and sinks used 
in the simulation. Here we did not attempt high preci- 
sion, but rather used simple local sources to do a quick 
check on the spectrum. Our results, shown in Fig. [7| 
agree well with experiment to within our statistical er- 
rors. 



VI. SIMULATION TECHNIQUES FOR 
HIGHLY-SMEARED OPERATORS 



It would be highly desirable to create new unquenched 
gluon configurations using the HISQ action in place of the 
ASQTAD action. The use of such configurations, with 
lattice spacings of order a = 0.1 fm, would significantly 
reduce any residual worries about errors due to taste- 
exchange interactions. 

Such simulations are complicated, however, by the 
heavily smeared, reunitarized links in the action. The 
additional smearing and reunitarization have no effect 
on the cost of the quark-matrix inversions required when 
updating gluon configurations, but they complicate cal- 
culation of the derivative D a (x, /i) S q of the quark ac- 
tion S q [U] with respect to an individual link opera- 
tor U/j,(x). Derivative D a (x,fj,) is defined by 

/(e i£ ° T X(z)) = f(U, l (x))+e a D a (x^)f(U^ + 0(e 2 ) 

(41) 

for any function / of U^(x). We developed and tested 
both an analytic and a stochastic version of the derivative 
for this action [2^. 

The analytic version employs a unitary projection 

V^(VV r yW (42) 

to reunitarize each smeared link V . The main obstacle is 
then computation of the gauge derivative of the inverse 
square root. Using the product-rule identity for the ma- 
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trix M = VV\ 

Mi (DAf-i\ + (dM~^ = -M~i (DM) M~i, 

(43) 

we can solve for D(VV 1 ' )~2 directly, both iteratively by 
the conjugate gradient algorithm, and exactly by first 
diagonalizing VVt 

Using the chain rule, we combine D(VV^)~^ with 
standard derivatives of the base action and of the 
smeared links V. We encode derivatives of the action and 
smeared links generically, allowing for run-time changes 
independently in either. The additional cost of comput- 
ing and combining D(VV^)~? with these is minimal. 

Several other analytic a ppr oaches have been developed 
for unitarized smearings |27l |2S[ . The authors of |27j 
addressed the problem of computing the derivative of the 
matrix inverse square root by replacing it with a rational 
approximation. In |28j |. the smearing itself is unitary, 
explicitly avoiding the need to reunitarize and the inverse 
square root. In both, the derivatives are then computed 
straightforwardly. 

In the stochastic approach, we define links 

U*(x) = e Ur >^U^x) (44) 

where r)(x,/j,) is now a field of traceless, hermitian 3x3 
random matrices, each with normalization 

(VijVlmtv = {^imSji - -5ij8 lri ^j , (45) 

and e is a very small number. Defining AS q (x,/i) to 
be the terms in S q [U v ] — S q [U] containing U^(x), the 
derivative can be computed efficiently from 

2eT a D a (x, n) S q w (r)(x, fi) AS q (x, fi)) n (46) 

averaged over a finite number of sets of random r/s. Our 
numerical experiments suggest that 10-100 sets of 77s are 
adequate for actions like HISQ. Wc will describe both our 
stochastic and analytic techniques in a later paper. 

VII. CONCLUSIONS 

In this paper, we have demonstrated that taste- 
exchange interactions are perturbative and we have 
shown how to use Symanzik improvement to create a 
new staggered-quark action (HISQ) that has greatly re- 
duced one-loop taste-exchange errors, no tree-level or- 
der a 2 errors, and no tree-level order (am) 4 errors to 
leading order in the quark's velocity v/c. The HISQ ac- 
tion addresses one of the fundamental issues surround- 
ing staggered-quark simulations by allowing us to esti- 
mate taste-exchange interactions through comparisons of 
HISQ with ASQTAD results. We presented numerical 
evidence that taste-exchange interactions in HISQ con- 
tribute less than 1% to light-quark quantities, like meson 



a = 1/8 fin a = l/llfm 
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FIG. 8: Masses for two different tastes of the D 3 and one 
for the D* using HISQ c quarks on lattices with lattice spac- 
ings a = 1/8 fm (e = -0.21) and 1/11 fm (e = -0.115). For 
HISQ, D 3 s are given (from left to right) for the 0-link and 
1-temporal-link mesons; only the 0-link splitting is given for 
the D 3 . The dashed lines indicate the results from experi- 
ment. The error bars shown are from the simulated values for 
bid, — m,, c /2 and — m r)c /2, since these quantities are 

insensitive to the c mass and the c mass was tuned to make 
the simulated value for m Vc exact. 



masses and decay constants, at lattice spacings as large 
as 0.1 fm. 

The suppression of all order (am) 4 errors by powers 
of (v/c) 2 makes HISQ the most accurate discretization 
of the quark action for simulating c quarks on current 
lattices. We demonstrated this with a new lattice QCD 
determination of the tp — rj c mass splitting, which agrees 
well with experiment. This result could be improved 
by computing the coefficient cms in the correction SCms 
(Eq. pg|)l to the HISQ action. 

Our final HISQ action is defined by Eqs. l|32H3B|) . We 
showed that the tree-level value (Eq. J29) for the Naik- 
term parameter e is adequate at the level of 1% errors for 
lattice spacings at least as large as 1 /8 fm. 

The HISQ formalism will be most useful for D physics 
and works very well for such mesons. In Fig. [8] for exam- 
ple, we show the D*—D s spin splitting from simulations 
with HISQ c quarks using both our 1/8 fm and 1/11 fm 
lattices. As in the charmonium case, theory and experi- 
ment agree well here to within errors. Another very sen- 
sitive test of our simulations is to compare computed and 
experimental values for the mass difference 2ttld s — m Vc , 
which is quite insensitive to tuning errors in the c mass. 
We obtain 956(14) MeV using HISQ on the 1/11 fm lat- 
tice, and 978(20) MeV on the 1/8 fm lattice. Both val- 
ues agree well with the splitting 956(1) MeV from exper- 
iment. Finally, we measured the speed of light and found 
c 2 = 1.00(4) for the D s on the coarse lattice (a = 1/8 fm 
with e = —0.21), confirming that the same coupling con- 
stant values work for both D s and r\ c mesons. It is impor- 
tant to appreciate that there are no free parameters avail- 
able for tuning in the extraction of any of these results; 
all QCD parameters were tuned using other quantities. 

As we will discuss in a later paper, the HISQ formalism 
is particularly useful for accurate calculations of quanti- 
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ties like fo, /d s5 D — ► 7rZf, and so on. These all require 
currents and these currents, even though they are con- 
served or partially conserved, have order a s (am) 2 renor- 
malizations. Consider, for example, cc annihilation into 
a (virtual) photon. The electromagnetic current can be 
computed by inserting photon link operators, 



C/Q ED ^exp 



X + fl 



dx ■ A^ ED j , 



(47) 



into the HISQ action and then expanding to first order 
in A QED to obtain 



• AQED tQED 

1 (U ° /J 



(48) 



Current jQ ED is not rcnormalized, because of QED's 
gauge symmetry, but it is not the only operator that con- 
tributes to cc annihilation. In addition there are 0(a 2 ) 
correction terms that contribute, 



ci a s a 2 m^a •eF QED ^ 

+ c 2 a s a 2 ^(C-eF) QED -7^, 



(49) 



just as for the gluon fields. The first correction term 
renormalizes the quark's magnetic moment, the second 
its charge radius. Such terms arc normally negligible 
since the extra derivatives introduce extra powers of (ap) 2 
and p is small. For c quark annihilations, however, the 
extra derivatives become c quark masses instead of small 
momenta, and so are more important. It is easy to com- 
pute coefficients c\ and c 2 in perturbation theory, but 
contributions to annihilation from these operators are in- 
distinguishable from those coming from the leading op- 
erator jQ ED . Consequently we can omit the corrections 
and, instead, introduce a renormalization constant for 

jQED. 



where 



jQED, continuum _ g ff jQED,HISQ ^q-j 



Z cS = l + ca s (am) 2 + 0(a 2 (am) 2 ). (51) 



The situation is very similar for (partially-conserved) 
weak-interaction currents. Again the current de- 
rived from the action requires 0(a s (am) 2 ) corrections. 
Once these one-loop corrections have been calculated 
and included, however, the remaining terms are only 
0(a 2 (am) 2 ). Consequently one- loop radiative correc- 
tions are all that is necessary to achieve 1-2% preci- 
sion for /d, D — > wlv, and similar quantities. Other 
discretization errors in these quantities should be of or- 
der 1% or less. With one-loop renormalizations, lattice 
results will be more accurate than current results from 
CLEO-c and the b factories. 
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APPENDIX A: GAMMA MATRICES 

The complete set of spinor matrices can be labeled by 
a four-component vector consisting of Os and Is (i.e., 

£ Z 2 ): 



7 „ = n w*- 



(Al) 



We will also sometimes use the more conventional, but 
equivalent (up to a phase) hermitian set: 

1, 75 = -jaxjylz, Jn, 

75^ = ^757^ 7^ = fbv,7*d 1 - ( A2 ) 
The 7„s have several useful properties: 

• orthonormal: 

Tr(7t 7m )=4(5„ m (A3) 

• closed under multiplication: 

7n7m = (-I)"'" 1 " 7 n+ m (A4) 

where 

= ^ m„ mod 2 (A5) 



v<p, 



and 



with 



n> = n v mod 2; 



(A6) 
(A7) 



v>p, 



(A8) 



• hermitian or antihermitian: 

In = In — In i 

• commuting or anticommuting: 

7n7m = (-1)^ 7 m In (A9) 

where 

m^ = m> + m< — m v mod 2 

jm u if m 2 even 



(nip + 1) mod 2 if m odd 



and 



m ■ n = n ■ m, 



(All) 



We thank the MILC collaboration for sharing their 
configurations with us. This work was funded by grants 



m = m (m ■ m) mod 2 = 0; (A12) 
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• permutation operator: if one uses the standard rep- 
resentation for gamma matrices, where 



7o = 



and 



1 

-1 




&3 = 



li = 



<?2 = I 

1 

-1 




(A13) 



-i 

1 



(A14) 



then there is at most one nonzero element in any 
row or any column of any 7„, and that element is 
±1 or ±i. Thus multiplying a spinor ij> by 7„ simply 
permutes the spinor components of tp, multiplying 
each by ±1 or ±i. 

A convenient notation, reminiscent of 75^, is 

7nm=7«7m- ( A 15) 

APPENDIX B: STAGGERED QUARKS 

The naive discretization of the quark action is formally 
equivalent to the staggered-quark discretization. Stag- 
gering is an important optimization in simulations; it is 
also a remarkable property. Consider the following local 
transformation of the naive-quark field: 

tl>(x)^Sl(x)x{x) ^{x)^x{x)tf{x) (Bl) 



where 



fi(x) = j x 



ft=0 



(B2) 



and we have set the lattice spacing a = 1 for convenience. 
(We will use lattice units, where a = 1, in this and all 
succeeding appendices.) Note that 



fi(x) = 7„ 



for rii. 



mod 2; 



(B3) 



there are only 16 different S7s. It is easy to show that 



a^(x) = f2 t (x)7 At r2(x ± p) 
1 = &(x)n(x) 



(-If,, (B4) 
(B5) 



xo + Xi + • • • + x M _i (see Appendix lA*j l. 
Therefore the naive-quark action can be rewritten 



where x5 

T L 



t/)(x)(7 • A + m)ip(x) = x(x)(a(x) • A + m)x(x). (B6) 

Remarkably the x action is diagonal in spinor space; each 
component of x is exactly equivalent to every other com- 
ponent. Consequently the x propagator is diagonal in 
spinor space in any background gauge field: 



(x(x)x(y)) x = s(x,y) l s 



(B7) 



where s(x,y) is the one-spinor-component staggered- 
quark propagator. Transforming back to the original 
naive-quark field we find that 



St- 



(V(x)v(y))v- = s(x,y)Q(x)^(y). 



(B8) 



This last result is a somewhat surprising consequence 
of the doubling symmetry. It says that the spinor struc- 
ture of the naive-quark propagator is completely inde- 
pendent of the gauge field. This is certainly not the case 
for individual tastes of naive quark, whose spins will flip 
back and forth as they scatter off fluctuations in the chro- 
momagnetic field, for example. The sixteen tastes of the 
naive quark field are packaged in such a way, however, 
that all gauge-field dependence vanishes in the spinor 
structure. 

Doubling symmetry is immediately evident in the stag- 
gered action, Eq. I|B6|) , since the action is invariant under 



X (x)^S f (0)x(x), 



(B9) 



which merely scrambles the (equivalent) spinor compo- 
nents of xi x )i and since 



B c (x)n( x ) = n(x)B c {o). 



(BIO) 



In simulations one generally discards all but one spinor 
component of x(x), resulting in highly efficient algo- 
rithms. The Symanzik improvements discussed above are 
trivially incorporated. 

Note, finally, that Eq. IjBlOfl implies 



B c (x)Cl(x)tf(y)Bl(y) = n(x)n+(y). (Bll) 

and therefore the naive-quark propagator, Sf, satisfies 

S F (x,y) = B c (x)S F (x,y)Bl(y) (B12) 

for any background gauge field. In momentum space this 
becomes 

S F (p,q) = B c (0)S F (p+(Tr,q+(Tr)Bl(0), (B13) 

which is an exact relationship that is useful in pertur- 
bative calculations. This last relation, which is easily 
checked at tree level (but true to all orders in a s ), shows 
that there is only one sixteenth as much information in 
the naive-quark propagator as naively expected. 

APPENDIX C: TASTE, NAIVE VS. STAGGERED 

The quark field ^(x) will typically contain contribu- 
tions from all 16 tastes. One can separate out (approxi- 
mately) the different tastes by blocking the field on hy- 
percubes that have two sites per side. One way to project 
out the C — taste, for example, is to average over the 
hypercube: 



(CI) 
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where Xb, with xb^ mod 2 = 0, identifies the hypercube, 
and the sum is over all 16 sites in the hypercube. Any 
component of ip(x) that has momentum p « C,tt with 
Q 7^ will be strongly suppressed by the average. (The 
suppression is 0(a 2 p 2 ) for a mode with momentum p + 
(it.) A £ ^ piece of ip(x) can be isolated by applying 
a doubling operator Bq to transform that component to 
C = 0, and then, again, averaging over the hypercube to 
isolate that component: 



1 £ 



— y, &d x B + Sx)ip(x B +Sx). 



The 16 blocked fields v4°, with one for each describe 
the 16 different tastes of quark, each now in the low mo- 
mentum sector (i.e., with p^ < 7r/2). 

With staggered quarks, one keeps only one component 
of x( x ) since the components are equivalent and decou- 
ple: 



(C2) 









(C3) 



We can use our hypercube blocking (Eq. l|C2|l ) to trans- 
late this single field back into blocked fields for "ordi- 
nary" quarks of different tastes £: 



Bq(xb + Sx) ip(xB + Sx) 



:^O(fe)S C (0) 



Xi{xb + Sx) 






(C4) 



Only four of the sixteen blocked fields are independent, 
however, because B^(0)x( x ) is always proportional to 
one of only four different spinors: 



(C5) 



(This is because B^(0) = 7^ and all 7„s when applied 
to a spinor merely permute the elements of the spinor, 
multiplying each by ±1 or ±i.) Consequently we can 
reconstruct just four tastes, t = 1 ... 4, of blocked quark 
from a single staggered field: 



Xi( x ) 






















Xi(x) 






















Xx(x) 






















. Xi{x) _ 



^r!(fe)x (t) xi(^ 

Sx 



Sx) 



(C6) 



where the x^'- 1 are un it spinors with xP = ^i,t- 

This last formula defines the standard taste basis for 
staggered quarks. It is of more use conceptually than 
practically, and we will not need it in what follows. 



Formula (|C4|I shows that some £'s become indistin- 
guishable in the staggered approximation. In fact, by 
explicit calculation (given our specific representation of 
the 7 matrices), one finds that the £s fall into four equiv- 
alence classes composed of indistinguishable £s: 



.4 
B 
C 
D 



0000 0001 0110 0111 

0010 0011 0100 0101 

1000 1001 1110 1111 

1010 1011 1100 1101 



(C7) 



These four classes correspond to the four tastes of stag- 
gered quark; all the Cs in a single class give the same 
staggered- quark field -ip^ (Eq. (|C6jl l from Eq. ljC4(l . 

The equivalence of the £s within a single class is pre- 
served under addition of £s: for example, adding any two 
Cs from class B always gives a(in class A, while adding 
any vector from class B to any vector from class C al- 
ways gives a vector in class D. The full addition table 
for these classes is: 



+ 


A B 


C D 


A 


A B 


C D 


B 


B A D C 


C 


C D 


A B 


D 


D C 


B A 



This kind of structure is necessary for the four-fold re- 
duction in the number of tastes due to staggering — that 
is, 16 tastes of naive quark are reduced to 4 tastes of stag- 
gered quark because we can consistently identify certain 
corners of the Brillouin zone with each other in the stag- 
gered case. 



APPENDIX D: NAIVE-QUARK CURRENTS 

Naive quarks lead to huge numbers of nearly equivalent 
mesons. Where in ordinary QCD one might consider a 
single meson operator, J„ = ijjj n ip, one has 16 point-split 
operators in the naive theory each of which couples to a 
different meson: 

4*\x) ^ ^{x)^H{x) 

= ip(x)^ n tp(x + dx sn ) 

« x( x )lsXi x + Sx sn ), (Dl) 

where s is any one of the 16 four- vectors consisting of 0s 
and Is only, 



Sx sn = (s + n) mod 2, 



(D2) 



and link operators are implicit (and a = 1 still). Each 
of these operators creates a different version of the J n me- 
son; they are orthogonal. This is because 



4 s > - (-i) s « 4 s ) 



(D3) 
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under a doubling transformation where 

il>(x) - B c=fl (x)^(x) $(x) -> ^(as)Bj(a:). (D4) 

Thus the four-vector s determines the bilinear's transfor- 
mation properties under arbitrary doubling transforma- 
tions; it specifies the bilinear's "signature" under dou- 
bling transformations. Since doubling transformations 
are symmetries of the naive theory, the doubling signa- 
ture is conserved: for example, 

(J«jWt) = ifr^s, (D5) 

which proves that each of our point splittings creates a 
different meson. 

Different signatures correspond to different variations 
of the same continuum meson, typically with slightly dif- 
ferent masses, etc. These are the different tastes of the 
meson. We label different tastes by the corresponding 
signature s in the meson's rest frame. 

Additional mesons are made by boosting particles into 
other corners of the Brillouin zone: 

Ptot ~ C 71 " (D6) 

for one of the sixteen £s. Such mesons would be highly 
relativistic in the continuum, but here they are equivalent 
to low-energy mesons because of the doubling symmetry: 
for example, the antiquark in the meson might carry a 
momentum near zero, while the quark carries C 71 " but is 
then equivalent to a zero-momentum quark state through 
the doubling symmetry. Pushing such momenta through 
a naive-quark bilinear, using for example 

X 

changes the quantum numbers of the meson created by 
the operator. To see this, consider current 

4^(X) EE {-lf^WfyPMx) 

= ^(x)Bl(x)^i/j(x) (D8) 
« x(x)j^"fsX(x + 6x sn ) (D9) 

which has signature (£ + s) mod 2. This current obvi- 
ously carries momentum p to t — C 71 " if averaged over all x. 
It is easy to show that flavor-nonsinglet mesons created 

(C s) 

by JA with different (s are all identical, and therefore 
all carry taste s. For a ud meson, for example, we can 
use the separate doubling symmetries of the u and d to 
prove that 

(Jp)(x)Jp)t(y)> = (uBl^dd^Bcu) 
= (VyUddrfPu) 
ee (JP(x)JPHv)) (D10) 

for all (s. 



The labeling on Jn is intuitive (and therefore use- 
ful) only if the current is averaged over x in such a way 
that ptot — C 71 " + P where —it/2 < p^ < 7r/2 for all /z; 
the combination of all £s and ps covers all of momentum 
space, so nothing is lost by this restriction and double 
counting is avoided. We can enforce this restriction on 
the momenta by redefining the current on a blocked lat- 
tice with one site at the center of every 2 4 hypercube 
on the original lattice, just as we did for the quark field 
(Section 0: 

J^ee 4 C ' s) (xb+6x) (Dll) 

where xb identifies the hypercube (with xb^ mod 2 = 
0). The blocked current creates a meson of taste s with 
momentum in the corner of the Brillouin zone. 
Note that momentum and taste conservation imply 

£<4 c '%)4 CV)t (y)) = o (D12) 

unless C — C an( i s — s'. Consequently the different 
operators are orthogonal when analyzed in momentum 
space. Typically we sum over space, however, but not 
time. In that case operators with the same signature, 
but where £' — C = ±i, can mix. This mixing leads to 
components in meson propagators that oscillate in time, 
as we discuss below. 

To summarize, there are 16 sets, each set labeled by 
taste s, of 16 identical mesons, each meson labeled by £, 
for each flavor-nonsinglet meson in the continuum. Each 
corner £ of the Brillouin zone has a single representative 
of each taste of meson. 

Flavor-singlet mesons are slightly more complicated. 
In a uu meson, for example, the quark and antiquark 
are created by the same field, and so do not have sepa- 
rate doubling transformations. Therefore the argument 

(C s) 

relating J„ s with different £s doesn't work. The only 
contributions that spoil this argument are from annihila- 
tion, where the meson's quark and antiquark annihilate 
into gluons. Annihilation gluons that contribute to, for 
example, 

£(Jp)(x)Jp)t (0)> (D13) 

x 

carry total momentum Qtt and so are far off-shell unless 
C = 0. This has two implications: First, annihilation 
contributions will be different for different (s. And, sec- 
ond, only the C = case h as the correct coupling between 
the flavor-singlet quarks and purely gluonic channels. In 
fact, only the taste-singlet state, among the £ = states, 
couples to the gluons since 

(4 s Hx))^ ex Tr( 7s ) 

= unless s = 0. (D14) 

A similar condition applies to nonzero £s as well. For 
each ( there is only one taste that can couple to glu- 
ons. The gluons are highly virtual if £ is nonzero. These 
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last contributions are taste- violating because they change 
quark taste along quark lines; they are removed by the 
contact terms discussed in this paper. 

It is not surprising that the flavor-singlet mesons are 
more complicated. They usually have to be, particularly 
in the pseudoscalar channel where the U(l) problem must 
be resolved. In our naive-quark theory, only the £ = 
s = mesons couple properly to gluons. The masses 
of pseudoscalars with ( = s = arc shifted properly by 
instantons in the chiral limit, so that the t/(l) problem is 
resolved. The £ = s = neutral pion is also the only pion 
that decays to two photons. The corresponding axial- 
vector current is only approximately conserved, even in 
the chiral limit, so anomalies are not needed to mediate 
the photon decay. 



APPENDIX E: STAGGERED-QUARK 
CURRENTS 



The 256 different mesons created by the naive-quark 
bilinears J Tl include 16 identical copies of each distinct 
taste of meson. Staggering the quark fields discards iden- 
tical copies, leaving just 16 distinct tastes. The naive- 
quark bilinears corresponding to the staggered-quark bi- 
linears are the ones that become diagonal after they are 
staggered. Any other operator creates mesons that mix 
different spinor components of the staggered quark opera- 
tor x{ x ) (Eq. HB1|) 1. and so is discarded when we stagger. 

We can identify the bilinears that survive staggering 
by noting that 



-> XX 



(El) 



after staggering provided £ = s. The staggered-quark 
operator is therefore 

jCM) = (-l)**tf(x)it>y n il>(x + 5x m ) 

= (3^{x)x{x) X {x + 8x sn ) (E2) 

where again Sx sn = (s + n) mod 2 (with a — 1), and 



fa\x) = \ Tr ( 7 hhn7: 



x-\-s-\-n ) 



(-1) 



x-{s < +n > ) 



n- (s+n) * 



(E3) 



Each taste s in the staggered theory corresponds to a 
specific corner of the Brillouin zone, withptot ~ stt, in the 
naive-quark theory. These operators have zero signature 
and so are unchanged under doubling transformations 
(Eq. ®). 

It is useful for formal analyses, though less so for sim- 
ulations, to introduce a slightly different definition for 
these bilinears by defining a new operator on our naive- 
quark field: 

7„ © £ s ip{x) = (-If 7] 7n tp(x © (n + «)), (E4) 



where © adds vector n + s to x "modulo" the hypercube 
that x lies in — that is, 

(x 9 n) p = x B u + (Bp - x Bp , + rift) mod 2 (E5) 

when x is in the hypercube labeled by site xb (with 
iBfi mod 2 = 0). Thus x © n is in the same hypercube 
as x. With this definition, we can redefine the staggered- 
quark current with spin n and taste s to be: 



ip(x)-y n (g>£ s il)(x). 



(E6) 



The © means that the new operators 7n © £ s have a 
simple algebra. If, for example, 

i>(x) = j n ^^t/j(x) (E7) 

then, using Eq. (IA9I) . 

7m © £r Ipfa) = 7m In © £r£,s ^(x), (E8) 

which implies that in general 

7m © £r 7n ® 6s = 7m7n®£r£s- (E9) 

Note that these definitions also imply that £s anticom- 
mute just as 7s: for example, 



(E10) 



We can effectively restrict our current to the C = s 
corner of the Brillouin zone by again blocking on 2 4 hy- 
percubes to obtain 

J s {n) (x B ) = — ^ XB + 5x ) 7« ® 6 i>{xB + Sx) 
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8x^ SZ2 

(Ell) 

where Xb^ mod 2 = 0. This formula is related to a more 
standard formula by staggering the meson operator: 



ip(x) 7 „ © 6 



X( x ) itll In lx+8x 3n X( x © Sx sn) 



(E12) 



where the shift 8x sn guarantees that the product of 7 s is 
proportional to the unit matrix. We can rewrite this in 
terms of a sum over Sxs in the positive unit hypercube: 

J2 X(x) 1 Tr ( 7 t 7 t 7n lx+Sx ) x (x © Sx). (E13) 
Averaging over the hypercube gives a standard formula 



for J t 



(«). 



^ x(xb+Sx)\Tt (jtjl x jnjsx>^ x(xb + Sx'). 

6x,8x' 

_ (E14) 
Taste and spinor structure in ip(x)"f n l S>^ s ip(x) are both 
described by the same kind of four vector consisting of 0s 
and Is. For this reason it is common practice to use the 
same terminology for describing taste as we do for spinor 
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structure. Thus, for example, %p(x)jz, ® £5 ^(a;) creates 
a pseudo-scalar meson with "pseudoscalar taste." In this 
case n = s = (1, 1, 1, 1), and, from Eq. (|E2|) . 

i>{x) 75 ® 6 ip(x) = (3 ( 5 ,) (x)x(x)x(x) 

= (-ir x x(x)x(x). (E15) 

A different pion is created by ip{x)^§ ^^{x), this one 
with axial-vector taste: 

^{x)^®^{x)=p£»\x)x{x)x{x®(i). (E16) 

The pion created by this operator is sometimes called a 
5 ® 5/i or "1-link pion" since the operator is split by one 
link or lattice spacing. Similarly 4>{x)^ <g>£,5ip(x) creates 
a 5 (g> 5 or 0-link pion. 

Other naive-quark operators can be recast in terms of 
the spinor/taste operators. For example, the naive-quark 
action has a chiral symmetry in the massless limit under 
transformations 



ip(x) -» exp(i6>75 <E> £5) i>(x) 
ip(x) -> ip(x) exp(i6»7 5 ®£ 5 ), 



(E17) 



where the £gs express the fact that the symmetry opera- 
tion does not translate (move) ip(x). 



Tree-level contributions come from Fig. but these 
vanish (by design) in the ASQTAD action because of 
ASQTAD's quark-gluon vertex. The only one-loop con- 
tributions that arc nonzero are those shown in Fig. [3 all 
other one-loop amplitudes have at least one quark line 
with only one gluon attached, and these vanish, again, 
because of ASQTAD's quark-gluon vertex. The inter- 
nal quarks and gluons in the remaining diagrams are all 
highly virtual, and therefore these contributions can be 
canceled by a sum of four-quark counterterms each con- 
sisting of a product of two quark bilinears (with one bi- 
linear per quark line). 

The types of quark bilinear that can arise from these 
diagrams are greatly restricted by color conservation, 
chiral symmetry, and the doubling symmetry of the la- 
grangian. Quark bilinears can only carry singlet or octet 
color: thus color structure is either ipip or ipT""^ where T a 
is an SU (3) generator in the fundamental representation. 
The standard chiral symmetry of the naive-quark action 
in the massless limit implies that only vector and axial- 
vector bilinears can arise in that limit: thus spinor struc- 
ture is either ip^j^ip or ^75 „ ip- Finally doubling symmetry 
(Eq. (|F1|) ~) requires that the bilinears in the counterterms 
must be point-split so that they are invariant under any 
doubling transformation of their quark fields (Eq. ©): 
that is, we have 



APPENDIX F: ONE LOOP TASTE-CHANGING 



ip(x)^ n tp(x + dx(n)) 



(F2) 



The Symanzik procedure for removing 0(a 2 ) one- loop 
taste-exchange effects from the ASQTAD action involves 
two steps: 1) we compute taste-changing amplitudes for 
qq — > qq with massless quarks in one-loop order us- 
ing lattice perturbation theory; and 2) we design local 
taste-changing counterterms for the staggered-quark ac- 
tion that cancel these one-loop amplitudes. We find it 
easiest to work with the naive-quark theory, converting 
to staggered quarks only at the end. 

In order a 2 we need only consider quarks at thresh- 
old — that is, quarks with momentum p = ±£71 for one 
of the sixteen £'s with £ M £ Z 2 (where, again, a = 1). 
Taste-changing amplitudes are ones where the incom- 
ing and outgoing momenta along any particular quark 
line differ by £-7r for one of the £s. For qq — > qq, over- 
all momentum conservation demands that the opposite 
change occur along the other quark line. Consequently 
while taste changes along individual quark lines, total 
taste is conserved. It is enough to compute amplitudes 
A(0, 0; £71", — Ctt) where the initial quarks have zero mo- 
mentum and the final quarks have momenta £7r for any £. 
(Recall that and — C 71 " are the same momentum on the 
lattice.) Amplitudes with other quark tastes in the ini- 
tial state are related to these amplitudes by applying the 
doubling symmetry to each quark line separately: 

A(Cw, C'V; (C + CK (C" - 00 = A(Q, 0; C?r, -Ctt) (Fl) 

where C;,C;'eZ 2 . 



tp{x)-j n T a ^{x + 5x{n)) 



(F3) 



where n is fi or 5/i, and Sx(n) is chosen to make the 
bilinear invariant under tp — > B^tp for all £. This last 
restriction implies that the bilinears must be staggered- 
quark operators Jn , as defined in Appendix El 

The one-loop amplitude, A(0, 0; C"", — (n), for a given £ 
is canceled by a specific set of counterterms, 

\ d « } l^« S) I 2 + color octet version ( F4 ) 



where s = £ (from Eq. l|E2|l '). Consider, for example, the 
case where = 1 for some fj, while £„ = for all v 7^ /i. 
The taste is s — ( or 5//. The spin n depends upon the 
spinor structure of A. For this £ there are four differ- 
ent spinor structures, each corresponding to a different 
staggered-quark operator: 



(VnvV') 2 
k^iA) 2 



. {J m )3 
1 -7(5^2 

. (Ji (5p) )2 



(F5) 



Each of these operators comes in color octet and color 
singlet versions, so eight counterterms are required to 
cancel A for £ 2 = 1. 
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The coefficients d„ ' in the counterterms are com- 
puted by on-shcll matching of the scattering ampli- 
tude A(0, 0; (tt, — (ir) to the sum of counterterms for 
each £. Our results are summarized in Table [H] Only 
the bilinears shown in Eq. I|28|) are needed here. The chi- 
ral perturbation theory devised by Lee and Sharpe pfll ] 
has ten additional bilinears, but eight of these are not 
invariant under doubling transformations for each quark 
line separately, and two do not involve taste exchange. 

APPENDIX G: MESON PROPAGATORS 

The fact that the single naive-quark field encodes six- 
teen different, identical tastes of quark has practical im- 
plications for simulations. For example, the operator 

J 5 (x) =^(x) l5 ^ b {x), (Gl) 

where Vff, is an unstaggered &-quark field and ip is a light- 
quark field, couples only to J p = 0~ mesons in the con- 
tinuum, including the B. If ip is a staggered field, how- 
ever, it also couples to + mesons [30|. 

The doubling-symmetry formula, Eq. ©, is useful in 
decoding such situations. The second contribution from 
J5 arises from high-energy states that couple to it. In 
simulations one normally forms correlators like 

G 55 (t) = £<Q|J B (x,t)4(0,0)|0), (G2) 

X 

where the sum over x guarantees that total three momen- 
tum p = 0. The operators, however, are not smeared in 
time, and so can create arbitrarily high-energy states. 



The 6-quark resists large energies, as these drive it far off 
shell, but the staggered quark is on-shell when its energy 
E w or when E w it. These two possible states of 
the light quark correspond to two different tastes. Con- 
sequently J5 couples to two different mesons: one whose 
light quark has taste £ = 0, and one whose light quark 
has taste C = (1, 0, 0, 0) or E w tt. 

The first of these two meson states is the normal 
0~ B meson. To interpret the second state, we trans- 
form the high-energy staggered-quark field back to a low- 
energy field (which we understand, since it behaves nor- 
mally) using the doubling symmetry formula, Eq. ijrjjl: 

^I^^OrX^ToX-l)'. (G3) 

Substituting in J5 we see that this component of the op- 
erator is 

^76*iU~» ~* ^(*757o)75*6 (-1)* = -#7o*b(-l)'- 

(G4) 

The operator tjri'yo^b, now with low-energy fields, couples 
to + mesons. (It is J = because there is no three- 
vector index. It is P = + because P70P = 70 where 
P = 70 is the parity operator.) The full correlator has 
two components: 

G 55 (i) -|(0|^ 75 * 6 |0-)| 2 e- s -* 

- (-I)*|(0|^ 70 vI/ b |0 + )| 2 e- £ + t . 

The second component, rather unconventionally, oscil- 
lates in sign from time step to time step. 
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